The effect of curvature and topology on membrane hydrodynamics 
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We study the mobility of extended objects (rods) on a spherical liquid-liquid interface to show 
how this quantity is modified in a striking manner by both the curvature and the topology of the 
interface. We present theoretical calculations and experimental measurements of the interfacial fluid 
velocity field around a moving rod bound to the crowded interface of a water-in-oil droplet. By using 
different droplet sizes, membrane viscosities, and rod lengths, we show that the viscosity mismatch 
between the interior and exterior fluids leads to a suppression of the fluid flow on small droplets 
that cannot be captured by the flat interface predictions. 



The dynamics of mobile inclusions in lipid membranes 
are fundamental to a variety of biological processes, in- 
cluding signal transduction [1 and the endocytosis of 
bacterial toxins [2 . Membrane inclusions, such as pro- 
teins or lipid "rafts" [3], are in many cases signifi- 
cantly larger than the lipids making up the membrane, 
so their dynamics can be studied in terms of macroscopic 
objects moving in a continuum fluid environment. Ad- 
ditionally, elucidating the mobilities and hydrodynamic 
interactions of colloidal particles at a fluid-fluid interface 
has important technological ramifications for the design 
and formation of membranes and capsules composed of 
particles assembled on droplets [5j. 

Low Reynolds number hydrodynamics in viscous mem- 
branes or interfaces differs substantially from the better- 
known problem of hydrodynamics in three dimensions. 
Because of the coupling between the two-dimensional 
membrane and its surrounding viscous solvent, in-plane 
momentum in the membrane is lost to the surrounding 
fiuid. The flows induced in the surrounding fluid gen- 
erate nonlocal couplings between the membrane velocity 
and stress. The net result of these effects is to introduce 
an inherent length scale - the Saffman-Delbriick (SD) 
length - which is set by the ratio of the (2D) membrane 
viscosity rym to the (3D) fluid viscosity t?, ~ '7m/'? ISJ. 
For cellular plasma membranes ~ l^m [7]. In con- 
trast, low- Reynolds number hydrodynamics in 3D fluids 
is a scale-invariant theory. The existence of an inherent 
length scale in membrane and interfacial hydrodynamics 
has complex and rather subtle effects on a variety of prob- 
lems, including protein diffusion in cell membranes ^[S], 
the flow of monolayers through channels [9] , the dynamics 
of monolayer domains [10 , microrheology of fluid-fluid 
interfaces pTj, and the mobilities of rigid and flexible 
extended objects in membranes [12]. 

In this letter we explore the effect of nontrivial interfa- 
cial geometry and topology on the hydrodynamics of vis- 
cous interfaces. We find two principal results with broad 
implications for particulate transport on curved or spher- 
ical interfaces. The first is that the compact topology of 
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FIG. 1: (color online) (a) Schematic illustration and calcu- 
lated membrane velocity fleld of a point particle of radius a 
(green disk) at the north pole subject to a force F, with a pin- 
ning force Fo at the south pole. Here the interior and exterior 
fluids are identical: 77+ = rj-, R/£o ~ 0.1 and R/a = 100. (b) 
Image from a connected rod of paramagnetic PMMA colloids 
at the interface of a water-in-hexadecane droplet decorated 
with microparticles. 

a spherical interface fundamentally alters the nature of 
the 2D interfacial velocity field. On a sphere, any vector 
field must include at least two vortices [13]. In contrast, 
there are no such singularities in the velocity field on 
an infinite, fiat interface. Secondly, the curvature of the 
interface introduces a new length scale - the radius of 
curvature, R. This geometric length scale competes with 
£0 in determining the hydrodynamics of particles embed- 
ded in the interface. Thus, geometry plays a role in par- 
ticulate transport on par with viscosity. These results 
have important biophysical implications, such as the re- 
tardation of the diffusive transport of membrane-bound 
proteins in highly curved regions of the membrane. 

To develop a theory of the hydrodynamics of curved 
surfaces we ignore inertial effects and impose force bal- 
ance at the interface. Specifically, we consider the re- 
sponse of an incompressible (i.e. constant area) spherical 
interface of radius R and its surrounding fiuids to a tan- 
gential point force F = Fy applied at the north pole 
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of the interface, as illustrated in Fig. [TJa). Because of 
the curvature of the interface, the in-plane force balance 
equation must be written in a manifestly covariant form: 
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where g^"' is the contravariant metric tensor and Da is 
the covariant derivative; the Greek indices run over the 
polar and azimuthal angles 6, (f>, respectively, f'^''* = 
FS{9)y /{2ttR'^) is the external force density applied at 
the north pole. The term in brackets in Eq. (jlj is the 
viscous force density resulting from gradients in the in- 
terfacial velocity field Va', the last two terms are the vis- 
cous stresses due to solvent inside (ct~) and outside (cr+) 
the spherical surface, a^j — ■q± [Aw^ + Djvf^^ — P±5ij, 
where P±, ?y±, and are the hydrostatic pressures, vis- 
cosities, and velocities, respectively, of the solvents inside 
(— ) and outside (+) the sphere. To determine the solvent 
stresses on the interface, we must solve the incompress- 
ible Stokes equation inside and outside the sphere using 
the "stick" boundary conditions v*|^^^ = v. 

It is convenient to decompose the dynamical system 
into normal modes consisting of the combined flows of 
the interface and the interior and exterior solvents. The 
deformations of a 2D interface can be decomposed into 
bending, compression, and shear modes. However, the in- 
compressibility of the interface prevents compression and, 
when combined with the incompressibility of the interior 
fluid, bending modes of the interface. The remaining 
shear modes can be written in terms of the (manifestly 
covariant) curl of a scalar field, Wq, — eapD^'^, where e^p 
is the alternating tensor. The combined interface and 
solvent system is diagonalizible in a basis of spherical 
harmonics |14l 115] . In the region exterior to the sphere, 
we retain only those terms that vanish at infinity. This 
eliminates the solution that corresponds to the uniform 
center of mass translation of the sphere with respect to 
the surrounding fiuid. 

By expanding the delta functions in f''''* in spheri- 
cal harmonics and applying the in-plane force balance 
condition Eq. ([l]), we determine the amplitude of each 
normal mode of the combined interface/solvent system 
generated by the external force. The interfacial velocity 
is then given by 



V • e =~Vsin (/) ^ — CSC OPl (cos 9) , 
1=1 
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V • = -l^cos 0^ - [cot ePl (cos e)+Pf (cos 6)] , (3) 
1=1 

where V ^ F/ (47r77„i), P™{x) is the l-th associated Leg- 
endre function, and 
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The upper limit on the sums is defined below. 

In Eq. Q we have defined two lengths in analogy to 
the SD length: £± = 'ri^/ri±. For the case of a flat inter- 
face between two differing fluids, these two lengths enter 
symmetrically, so that the only one length scale controls 
the interfacial hydrodynamics: £q — rj-^J {rj- + 77+). Here 
the two lengths enter independently. 

The most striking manifestation of the effect of the 
interior/exterior solvent asymmetry occurs in the limit of 
a large interior viscosity, ry+ 3> ?y_ . This causes the I = 1 
term in Eq. ^ to dominate the sums in Eqs. ([2]) and (|3|. 
As a result, the external force at the north pole causes 
a rigid rotation of the interface and interior fluid. The 
opposite limit 77-)_ 3> rj^ will not have an analogous effect. 
However, for a small enough sphere, i? ^ the same 
rigid body rotation is observed. Thus, geometry alone 
can have a dramatic effect on interfacial hydrodynamics. 

In order to prevent the rigid rotation of the sphere, we 
apply a pinning force Fq at the south pole that forces the 
total fluid velocity to vanish at the point (see Fig. [1|l). 
Because of the linearity of the Stokes equation, the total 
response of the fluids and interfaces is simply the sum of 
the individual responses to each force. A typical solution 
for the interfacial velocity fleld on the sphere is shown 
m Fig. [l}i. The appearance of a vortex in the upper 
hemisphere is required by topological constraints; there 
is a similar one placed symmetrically on the back side of 
the sphere (not shown). 

The particle's mobility /i, defined by v = /iF, is given 
by [see Eqs. ^ and ^] 
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The finite particle radius a sets the upper limit l^ax — 
8R/9a of the sums in Eqs. ([2]), (|3|, and ([5]). It is deter- 
mined by the requirement that the Stokes mobility for a 
sphere of radius a is recovered for the case 77+ ~ rj^ with 
a vanishing interfacial viscosity (77m ~* 0) [11 . 

In Fig. [2] we plot the dimensionless mobility r/m/i for 
a particle at the north pole of a pinned spherical mem- 
brane as a function of R for a variety of interior vis- 
cosities r/_. In all cases, the flat-interface SD result 
JJmMsD ~ ln(^o/a)/47r (horizontal dashed lines) is recov- 
ered in the limit R 00. However, when R/£+ <C 1, 
VmlJ' ln{R/a)/2n (dot-dashed hue). Hence, particle 
mobilities in high-curvature membranes, R/£o ^ 1; are 
depressed relative to the SD result (this regime is not 
shown for the green/gray curve in Fig. [2]). For intermedi- 
ate curvatures, the mobility on a sphere can be enhanced 
or suppressed relative to the SD result, depending on the 
value of the ratio rj^/rj^. For rj^/rj^ < 1 (green/gray 
curve), the mobility in a spherical membrane is larger 
than the flat membrane mobility because the more vis- 
cous fluid is bounded inside the spherical membrane; con- 
sequently, it dissipates less energy than in the case where 
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FIG. 2: (color online) Dimensionless mobility rymM for a parti- 
cle at the north pole of a pinned spherical membrane as a func- 
tion of the membrane radius R, for 1+ — 10fj,m, a = O.Ol/xm, 
and 77- = 107^+ (green/gray curve), rj- = r]+ (black curve), or 
J7_ = 0.1?;+ (dotted curve). The dashed and dot-dashed lines 
indicate the theoretical asymptotic results (see text). 

its domain is unbounded. Conversely, when rj^/rj^ > 1 
(dotted curve), the mobility in a spherical membrane is 
suppressed relative to the flat membrane case. 

To calculate the mobility of extended objects, we use 
the superposition principle applicable to this linear sys- 
tem [16 . Specifically, we consider a rod of length L 
embedded in the membrane. We approximate the rod 
by iV -t- 1 disks of radius a separated by a distance b, 
where L = Nb + 2a. We also apply a pinning force 
at the south pole that sets the fluid velocity to zero 
at that point; this force mimics the adsorption of the 
droplet on the substrate in our experiments [l? . Using 
the superposition principle, the total interfacial velocity 

v'°\0, 0) = ES' Xc.M^, <l>^)■ Here, F'^) is the 
force applied to the disk at the point {6i,(j)i); i = cor- 
responds to the south pole, and i = 1, .., N + 1 labels the 
disks in the rod. We consider only forces parallel to the 
rod axis. The response function Xa./3(x;y) gives w)j°'(x) 
due to a unit force in the /? direction applied at y. 

To determine the forces F^*^ we require that the total 
fluid velocity vanishes at the south pole, and that each 
disk in the rod move with unit velocity. These constraints 
provide a set of -I- 2 linear equations that determine 
F*^*^. Summing the -I- 1 forces acting on the rod mov- 
ing at unit velocity gives the inverse mobility of the rod. 
Using this same set of forces we can also calculate the 
entire velocity field in response to the rod's motion, both 
on the sphere and in the surrounding fluids. 

We performed experiments to measure the flow flelds 
on spherical droplets coated with a monolayer of small 
(370 nm) sterically-stabilized polymethylmethacralate 
(PMMA) particles [IH]. Water droplets (77. = IQ-^ 
N s/m^), typically 30-lOO^m in diameter, suspended in 
hexadecane (77+ — 3.34 x 10"'^ N s/m^) provided the 
spherical interface. The PMMA particles served dual 
roles: to set the interfacial viscosity and to allow the 
flow field to be measured using video microscopy and 
particle-tracking software. To create the flow fleld, we 



added micron-sized paramagnetic polystyrene particles 
(carboxylate-functionalized, DVB-crosslinked, 0.95/i-m in 
diameter; item #MC04N, lot 3251 from Bangs Labora- 
tories). These paramagnetic spheres also adsorbed at 
the interface. In the presence of a magnetic field, they 
formed a single rod-like aggregate on the droplet. This 
rod was then moved along the surface of the droplet by 
a permanent magnet brought close to the sample. 

Samples were observed under bright-field microscopy 
using a Zeiss Axi overt 200. Images were captured at 30 
frames/second and analyzed with particle tracking code 
written in IDL [19]. The PMMA beads were first tracked 
without any magnetic field gradient present, and their 
mean squared displacements were used to determine the 
interfacial viscosity 77m. To measure the flow field, « 10^ 
PMMA particles were tracked while a single chain on the 
interface was moved at speeds of approximately a few 
fiva/s. For each droplet, the process was repeated twelve 
times, and the mean and statistical uncertainty of the 
flow velocities was measured. Droplets of different radii, 
chain lengths, and viscosities were used. 

In Fig. [3] we plot (points) the measured interfacial ve- 
locity, Vexp — Voxp{x)'y, along the line that perpendicu- 
larly bisects the rod (i.e. the line = 0,7r, < 6* < tt), 
as a function of the distance x = R sin 9 from the north 
pole in the x-direction. We also show the predictions 
of the flat interface theory (FIT, dashed lines) [12] and 
spherical interface theory (SIT, solid lines); each curve is 
obtained by direct calculation using no adjustable param- 
eters. We account for the rod thickness in the theory by 
setting Uoxp = vq everywhere within the rod. In Figs. [3]^ a) 
- Fig. [3f^c) , we show a sequences of droplets demonstrat- 
ing the increasing effect of curvature. In Fig.[3|a), where 
R '3> io, we see that both the FIT and SIT agree with 
the data. In Fig. [3|b), where R ^ io, but is now com- 
parable to L, the effects of curvature begin to be seen. 
However, only when £0 approaches R, as in Fig. ^c), 
does the effect of curvature become dramatic. Here the 
velocity field decays significantly more rapidly than FIT 
predictions. This is primarily because of the viscosity 
mismatch between the less viscous interior (water) and 
more exterior (oil) fluids, which leads to a relative de- 
crease in the particle mobilities on the spherical interface 
(dotted curve in Fig. |2|. In this case the SIT prediction 
is a significant improvement over that of the FIT. 

This work demonstrates the considerable effect of in- 
terfacial curvature and topology on the transport of par- 
ticles embedded in the interface. The topology of the 
sphere requires the formation of vortices in steady-state, 
zero Reynolds number flow, and the ratio of the radius 
of curvature of the interface to the SD length charac- 
terizing its 2D hydrodynamics determines the location 
of these vortices relative to the moving particle. These 
effects have direct implications for the kinetics of par- 
ticulate aggregates on the surface of droplets. In the 
future we will consider related problems on cylindrical 
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FIG. 3: Experimental data (points) and theoretical predic- 
tions (solid curves) for the interfacial velocity ?;cxp measured 
along the perpendicular bisector of the rod and normalized 
by the velocity vo of the rod. The dashed line shows the 
corresponding predictions for rod's velocity field on a flat 
interface [12]. The theoretical curves were generated using 
the parameters: (a) £o = 0.2fim, L — 6.4/^m, R = 17.5/im, 
A'' = 30. This data corresponds to Fig 1(b); (b) (.q = 0.3/im, 
L — 31.6/im, R — 54^m, A'^ — 35; (c) £o = 9.5/im, L = 7/im, 
R = 26^m, iV = 10. 



membranes, where we can separate the effects of curva- 
ture and topology. We have also demonstrated that a 
viscosity mismatch between the interior and exterior sol- 
vents can lead to either significant increases or decreases 
in the diffusivity of interface-bound particles depending 
on sphere radius. These effects may play an important 
role in protein transport on membranes separating the 
viscous cytosol 20J from extracellular fluids. 
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